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Introduction 


In the past two and one-half years of the NAGW-9 grant period (from 
November 1, 1983 to the present), we have rigorously pursued a research 
program in the theoretical and numerical modeling of solar activity and 
its effects on the solar atmosphere within the context of magnetohydro- 
dynamics. Specifically, our scientific objectives have been concerned 
with the physical mechanisms for the flare energy build-up and subsequent 
release. In addition, transport of this energy to the corona and solar 
wind was also investigated. The essence of this research program is to 
develop well-posed, physically self-consistent, numerical simulation 
models that are based upon magnetohydrodynamic theory. This approach 
enables us to understand the physics of flare energy build-up and its 
release as well as its effects on the coronal and heliospheric medium. 

Our approach is to conduct a systematic investigation of the basic prin- 
cipal processes that determine the macroscopic dynamic behavior of solar 
and heliospheric phenomena. In summary, we may claim that our scientific 
approach has met with some success --a total of twenty- three articles 
has been accepted and published in major refereed journals; Astrophysical 
Journal, Solar Physics, Astrophysics and Astronomy, Astrophysics and 
Space Sciences, Journal of Geophysical Research and also in some impor- 
tant conference and workshop proceedings. A detailed list is included in 
the Appendix. The major achievements are summarized according to sub- 
jects as follows. 

Evolution of Active Region 

A recent study by a number of investigators (Krall, et al . , 1982; 
Hagyard et al. , 1984; Hagyard, 1984a, b; 1986) has revealed that photo- 
spheric shear motion is one of the major causes of the flare. Based on 
this scenario, a two-dimensional, time -dependent, non-planar magneto- 
hydrodynamic (MHD) model in cartesian coordinates developed by this P.I. 
(Wu et al. , 1983a) was used for analyzing the evolution of active re- 
gions . Interesting findings are briefly summarized in the following: 

Evolution of Currents and Fields in an Active Region 

In this specific subject, we prescribe a shear motion at the lower 
boundary (i.e., photosphere), and then examine the atmospheric responses 
through the magnetohydrodynamics (MHD) process. Figure 1 shows; (a) the 
magnetic energy intensity; (b) magnetic field lines,; (c) parallel 
current (J ) intensity; and (d) perpendicular current intensity 1^200 
seconds after introduction of a photospheric shear motion (1 km s ) . 

The initial, characteristic plasma beta is taken as 3 “0.1, in which 
3 = lihl- D oh T 0 anc j corresponds to B q - 1500 gauss. From |hese results, we 

note that°the concentrations of magnetic energy (erg km ) and of elec- 
tric current (A km ) develop in similar ways. The location of these 
concentrations is in the neighborhood of the neutral line and in the 
lower parts of the atmosphere. Furthermore, the hill-valley type struc- 
tures of the magnetic energy intensity and the current intensity are also 
observed during this development stage (details, see Wu et al. , 1985). 

In general, the development of the magnetic energy and electric current 
intensity is much like that of twisted ropes. It has been shown theore- 
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tically (Sokolov and^Kosovichev, 1978) that magnetic field gradients of 
the order of 1 G km will trigger instability. Hence, we propose that 
the type of shear motion under study will lead to a catastrophic state of 
energy release which may be identified as a flare. In addition, the 
present self-consistent MHD model gives the physical parameters as func- 
tions of time and space and therefore can be tested by observation. 

Figure 2 shows the comparison between the build-up of excess magnetic 
energy in an initially pre-twisted (i.e., force-free) field and the 
build-up in an untwisted (i.e., potential) field configuration due to 
photospheric shearing motions. The excess magnetic energy build-up in an 
initially pre-twisted field is almost two orders of magnitude higher than 
in the case of an initially un-twisted field. 

Filament Formation Due to Photospheric Shear 

Another interesting physical feature has also been revealed by this 
numerical model (Wu^et al. , 1983a). Again the initial state is cljigsen ^n 
isothermal (T - 10 K) model in hydrostatic equilibrium (n - 10 cm ) ; 
and permeated by a dipole potential magnetic field that has a strength 
of 7 Gauss at the photospheric level (i.e., x-z plane). With this 
initial equilibrium state, we introduce a spatially-distributed sine 
curve type of shear motion (i.e., in the x-z plane) . The shear speed is 
zero at the neutral line and is a maximum (1.0 km s ) at the edge of the 
magnetic arcade. The evolution of the magnetic field configuration due 
to this prescribed shear motion is shown in Figure 3. The corresponding 
induced current distribution is shown in Figure 4a, b for several times at 
two specific heights (y - 1,500, 2,500 km) and distributed over a rela- 
tively small (8,000 km) scale in the horizontal direction (x-axis). 

Also, the induced current distribution as a function of altitude at a 
horizontal position of 2,500 km from the neutral line is shown in Figure 
4c. The corresponding plasma properties (i.e. density and temperature) 
at a specific height and horizontal position are given in Figure 5. To 
summarize these results, we make the following observations; 

(a) Figure 3 shows that the field lines have orderly changes 
up to 2,000 s after introduction of the shear motion in 
the x-z plane. A drastic deformation of magnetic field 
topology occurs during the interval between 2,000 s and 

2.500 s. 

Physically, this result shows that some of the field 
lines collapsed to form a current sheet (Figure 3) . The 
appearance of a high density region is coincident with 
the current sheet; it further shows that the outer edge 
of the magnetic field has locally been pushed open. 

(b) A further, specific, comment may be made about the 
current sheet. From Figure 4, we can identify the 
location at t^- 2,500 s, of the maximum current density 
(-10 Amp km ). Specifically, this peak is located at 

3.500 km from the neutral line and within the height 
interval between 1,500 km and 2,500 km. 
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(x = 2,500 km) with several times 
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Fig. 5. The pies., density (ale) and temperature (e,£,g) distribution as a 
function of height and horizontal axis at several times (1 000 s 
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(c) The density and temperature plots in this region (Figure 
5) show that the density has been enhanced 30% - 40% and 
temperature has decreased -20%. Based on these physical 
properties, we suggest that the shear motion induced a 
region which is characterized by a high current plasma 
density, high density and low temperature. This combina- 
tion of physical parameters contains all the character- 
istics for a filament. Since the MHD model for this 
simulation incorporates ideal MHD theory, we may con- 
clude that the filament could be formed by a plasma 
pinch effect. Details of this simulation have been 
given by Wu (1986) . 

Three-dimensional Non-linear Force-free Magnetic Field Model 

It has been recognized from observations that use of potential and 
linear force-free magnetic field models to interpret data for active 
regions is inadequate. A compromise approach - constructing a "patchwork 
quilt" representation of the field of an active region by combining 
magnetic fields derived from solutions of the linear equations for dif- 
ferent value of alpha - is of questionable value because it has no mathe- 
matical basis. It is certainly inappropriate for describing the evolu- 
tion of magnetic fields when important non-linear physical processes such 
as energy storage and release and MHD instabilities are involved. There- 
fore , we have developed a numerical model which be used for extrapolating 
non-linear force-free magnetic fields from a source surface, i.e., from 
observed vector magnetic fields at the photospheric level. A paper 
concerning the details of this numerical model was presented by Wu et al . 
(1985). Figure 6 shows the numerical results deduced from the observa- 
tions obtained by NASA/MSFC solar vector magnetograph (September 15, 
1982). Figure 6a shows the line-of-sight magnetic field (B ) contours 
and transverse magnetic field (B and B ) in vectoral form interpreted by 
using our non-linear force-free model a? the photospheric level (z - 0) , 
Figures 6b, 6c, and 6d show both the line-of-sight and transverse magne- 
tic field at various heights (-2,000 km, 10,000 km, and 15,000 km, 
respectively) . 

2.2. Coronal Dynamics 

Using the ideal MHD model with pressure pulses to study coronal 
dynamics has been done extensively. It was found that the current model 
could not give satisfactory physical interpretations of the observations. 
The main difficiency of this model is possibly due to an unrealistic 
initial state in the ideal MHD model with the pressure pulse. Further- 
more, the pressure pulse may not be the primary driving force. There- 
fore, we have begun to explore other physcial possibilities. Currently, 
we have initiated two new approaches and some preliminary results were 
obtained which we shall describe in the following: 
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Fig. 6. Extrapolated solar magnetic field based on nonlinear force-free model: 

(a) the observed magnetic field by MSFC vector solar magnetograph, 

(b) extrapolated nonlinear force-free magnetic field at 2,000 km height 

(c) extrapolated nonlinear force-free field at 10,000 km height and 

(d) extrapolated nonlinear force-free field at 15,000 km height. 





Magnetohydrostatic Equilibrium Atmosphere 


As we have pointed out earlier, the initial state of an isothermal 
hydrostatic equilibrium atmosphere with a potential magnetic field is 
inadequate, thus, we have constructed a magnetohydrostatic equilibrium 
atmosphere. In this sense, the initial state is now in a coupled state 
(i.e., non force-free state c.f., Wolfson, 1985) such that there is 
current flowing through the atmosphere. To obtain this solution, we have 
solved the complete time- independent MHD equation in two -dimens ions with 
the current structure as shown in Figure 7. Now we are ready to seek 
transient solutions with this newly established initial state which 
should enable us to study the dynamical evolution in the corona with a 
coronal streamer in magnetohydrostatic equilibrium. In the context of 
MHD theory, this initial state is very much different from the previous 
models (Nakagawa, et al. ; 1978, Wu et al. , 1978; and Steinolfson et al . , 
1978) . The present model has an initial state that couples the plasma 
properties and the magnetic field; the earlier models possessed an 
initial state that decoupled the plasma and the magnetic field 
properties . 

Formation of Plasmoid 

In this study, we use our ideal MHD model (Wu et al. , 1983), but the 
initial magnetic field topology is now changed to a helmet-steamer- type 
configuration as shown in Figure 8b. One can then postulate a lateral 
motion by applying a non-unifojm movement (with a maximum unrealistic 
velocity magnitude of 40 km s ) of magnetic field lines at the coronal 
base according to the pattern shown in Figure 8a. The reason we have 
chosen the unrealistic high velocity is that we chose an unrealistically - 
large physical domain (i.e. a hemisphere of the sun) for the purpose of 
this illustration. The initial plasma beta (g ) at the base is 0.1, 

L Q 

thereby implying strong magnetic control. Evolution of the magnetic 
field is shown in Figures 8c through 8f at times of 500 s, 1,000 s, 1,500 s 
and 2,000 s, respectively. Several interesting results can be noted 
and are summarized in the following: 

(a) The diverging motion from the central axis (i.e., the 
neutral line) leads to rising field lines in the closed 
field topology. This phenomenon has also been shown 
analytically by Low (1981) and numerically by Wu et al . 

(1986a). 

(b) As shown in Figure 8a, we also simultaneously introduce 
converging motion toward the axis of the open field 
configuration. As a result of this additional pertur- 
bation, high pressure regions develop on both sides of 
the neutral line. This high pressure acts as an 
"external" force that pushes the field lines together 
at a height greater than the neutral line. It appears, 
therefore, that there is a tendency for some of the 
field lines to be pinched off such that a plasmoid 
could be formed and subsequently, propagated outward. 

Actually, reconnection in this model is impossible 
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Fig. 7. A computed current distribution for a magnetohydrostatic corona. 





because the currently-considered model is an ideal MHD 
model. However, these results indicate that there is 
still another mechanism that could alter the field 
topology, even in the absence of a dissipative recon- 
nection mechanism. It is therefore, indicated that 
this additional consideration should seriously be taken 
into account when we interpret observed phenomena. 

Heliospheric Dynamics 

In this subject, our emphasis is centered on the propagation of solar 
originated disturbances into heliospheric space. One of the fundamental 
physical mechanisms in this context is the shock physics, i.e., struc- 
ture, acceleration and deceleration of shocks. In the following, we 
shall describe some interesting findings during this grant period. 

Shock Formation in the Heliospheric Space 

In order to illustrate the formation of interplanetary shocks (for 
details, see Wu, 1984), we have used the model developed by Wu et al. 
(1983b) with a rather weak disturbance (i.e. initial shock velocity being 
1.000 km sec -1 ) with a "piston” driving duration of 3,500 s. The three- 
dimensional representation of the disturbed solar wind configuration at 
various times is shown in Figure 9. Because the maximum strength of the 
disturbance (introduced at <J> - 30 ) is supersonic, the fast MHD shock 
propagates radially outward with heliolongitudinal spreading. At 
t - 20 hr, the early phase of heliolongitudinal spreading of the fast 
forward shock is beginning to appear. Formation of reverse MHD shock 
(caused by the finite pulse duration) is about to be seen from our choice 
of viewing angle in Figure 9a. At t - 40 hr, the fast reverse MHD shock 
begins to appear in addition to the heliolongitudinally- spreading fast 
forward MHD shock. The spreading phenomenon caused by non-linear wave 
coupling becomes more apparent as time progresses as shown in Figure 9c 
and 9d. 

Three-dimensional , Time -dependent Ideal MHD Model 


From our two-dimensional, non-planar, time -dependent MHD study (some- 
times referred to as 2-1/2 /d) , we have been convinced that a three- 
dimensional model could clear up a lot of ambiguity for the geometry 
which occurs in the two-dimensional non-planar case. Also, the great 
improvements of computer facilities has opened the possibility of real 
progress in this area. Thus we have started to develop a three- 
dimensional, time -dependent, ideal MHD model. We have met with some 
success, some of the preliminary results were published (Han et al . , 

1984); Dryer et al. ; 1986, Wu et al. , 1986b). Figure 10 shows represen- 
tative results obtained from our three-dimensional MHD model. This 
figure exhibits the magnetic field topology projected on the meridional 
(r - 0) plane and ecliptic plane (r - <j>) from 18 R (solar radii) to 76 R 
at various times (i.e., 5, 10, and 20 hrs. after initiation of the s 

perturbation). From these results, we observe a three-dimensional magne- 
tized plasma blob propagating outward due to the shock that is assumed to 
have been initiated from a flare. Also, we observe that the field lines 
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are pinched together in front of the blob, thereby creating a current. 
Since this is a full three-dimensional model, we would be able to esti- 
mate the volume of this blob of magnetized plasma on the basis of these 
results. The advantage of using the three-dimensional configuration is 
immediately apparent. 
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Magnetohydrodynamics (MHD) Modelling of Flare Energy Buildup, The Energy 
Release Phase, and Its Propagation into Heliospheric Space, by S. T. Wu and 
S. Panitchob in Proceedings of the SMA 1986. 

Magnetohydrodynamic (MHD) Modelling of Solar Active Phenomena Via Numerical 
Methods, by S. T. Wu, Solar Physics, 1986 (in press) 

MHD Simulation of Disturbed Solar Wind in Three Dimensions by S. T. Wu, M. 
Dryer and S. M. Han, Proceedings of the Symposium on Heliospheric Study in 
Japan (T. Watanabe, ed) December 1985. 

Filament Formation Due to Photospheric Shear by S. T. Wu and Y. C. Xiao, 
Proceedings of the CCP Workshop (A. Poland, Ed.) NASA SP Publication, 1986 
(in press) . 

Formation of Active Region and Quiescent Prominence Magnetic Field 
Configurations by C-H. An, J. J. Bao, and S. T. Wu, Proceedings of the CCP 
Workshop (A. Poland, Ed.) NASA SP Publication, 1986 (in press). 

A Linear MHD Instability Analysis of Solar Mass Ejections with Gravitation by 
M. T. Song, S. T. Wu and M. Dryer, 1986 (submitted to Solar Physics). 

A Three-Dimensional, Time -Dependent Numerical Modeling of Super-sonic, Super- 
Alfvenic MHD Flow by S . M. Han, S. T. Wu, and M. Dryer, July 1986 (submitted 
to Computers and Fluids). 
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Numerical Simulation of Flare Energy Buildup and Release Via Joule 
Dissipation by S. T. Wu, J. J. Bao, and J. F. Wang, J. of Adv. Space Research 
1986 (in press) . 

Numerical Tests of a Modified Full-Implicit-Continuous -Eulerian (FICE) Scheme 
with Projected Normal Characteristic Boundary Conditions for MHD Flows by S. 
T. Wu and J. F. Wang, International J. of Engineering Science, 1986 (in 
press) . 

Numerical Simulation of Mass Injection for the Formation of Prominence 
Magnetic Field Configurations by C-H. An, J. J. Bao and S. T. Wu, Solar 
Physics, October 1986 (submitted). 

The Use of An Adaptive Grid Finite Difference Method for Spherically 
Symmetric One -Dimensional Magnetohydrodynamic (MHD) Wave Propagation by S. 
Panitchob and S. T. Wu, January 1987, International J. for Methods in Fluids, 
(submitted) . 

The Observable Characteristics of Flare Energy Release II. The 
Interpretation of High Speed Soft X-ray Fronts by M. E. Machado, Y. C. Xiao 
and S. T. Wu, Ap.J.,1986, (submitted). 
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Travel 


The following is a list of travel by the P. I. and support personnel 
supported by NASA directly through this grant. 

1/83, S. T. Wu, to attend the SMM Workshop in Greenbelt MD and serve as team 
leader for the group studying flare energetics 

1/83, H. M. Chang, to attend the SMM Workshop held at Goddard Space Flight 
Center, Greenbelt, MD. 

5/83, S. T. Wu, to visit U. S. Air Force Space Divison and make a 
presentation on Solar Flares and its propagation into Heliospheric Model and 
Geomagnetic Storm Predictions for Operational Uses. 

6/83, James E. Dowdy, Jr., to attend and participate in the Solar Physics 
Division Meeting of the American Astronomical Society. 

6/83 S. T. Wu, to attend NASA/Solar-Terrestrial Physics Workshop as a 
panelist and NASA/SMM Workshop as Group Leader. 

8/83, S. T. Wu, to attend and present paper MHD Analysis of the Evolution of 
Solar Magnetic Fields in an Active Region" at the International Astronomical 
Union Syumposium #107. 

11/83 S. T. Wu, to attend International Solar Workshop in China as an advisor 
to the Chinese Academy of Science and present paper and to conduct a Solar- 
Interplanetary Seminar at Nat'l Centeral Univeriaty in Chengli, Taipei and 
Visit JPL and visit Stanford Univisty and Visit Lockheed Missile Research 
Laboratory in California. 

1/84 S. T. Wu, To attend AIAA 22nd Aerospace Science Meeting and present 
paper in Reno, Nevada. 

2/84, H. M. Chang, to attend 3rd meeting of the SMM workshop in College Park, 
Maryland. 

4/84 S. T. Wu, to attend and present paper at 6th annual technical conference 
ASME Solar Energy Division 

5/84, S. T. Wu, to attend and participate in 12th SE Conference on 
Theoretical and Applied Mechanics SECTAM XII. 

5/84, S. Panitchob, to attendand participate in SECTAM XII 


22 


6/84 S. T. Wu, To attend and participate in the Solar Terrestrial Prediction 
Workshop in Paris, France & present paper on Shearing Motions Within a solar 
Active Region with MHD Development of Magnetic Energy Buildup & Potential 
Flare Site Prediction. To attend COSPAR Planery Meeting/SMA Sympsoium in 
Graz, Austria representing the US as a memeber and session chairman of the 
SMA/COSTEP Steering Committee and present a paper on Solar Induced Geometric 
Disturbance. Visit Max Planck in Munich, Germany and Max Planck in Lindau, 
German, The Zurich Inst, of Techn. in Zurick, Sweden to present lectures and 
for consultation on Solar Physics Reserach and Von Karman Inst, of 
Fluliddynamics in Brussels, Belgium for consultation on Numerical methods of 
solar MHD problems . 

8/84, S. T. Wu, to consult with Dr. B. T. Tsurutani at JPL in Los Angeles, 

CA. 

9/84, S. T. Wu, to visit Goddard Space Flight Center to write SMM worksho 
report in Greenbelt, Maryland. 

9/84, S. T. Wu, to visit Space Enviornment Lab/NOAA for consultation in 
Boulder, Co and make presentation at AF Space Command in Colorado Springs, 

Co . 

11/84, S. T. Wu, to attend, participate in the U.S. -Japan Seminar in 
Heliomagnetosphere as chairman in Kyoto and letucurte on the subject of MHD 
instability at National Central Univerity in Taipai, and visit U. of 
California at Berkeley in San Francisco, CA. 

1/85, S. T. Wu to attend AIAA 23rd Aerospace Science Meeting and serve on 
Plasmadynamics and Lasr Technical Committee in Reno, Nevada 

3/85, S. T. Wu to visit AF System Command to make presentation and visit HAO 
for consultation. 

4/85, S. T. Wu, to attend CPP meeting in Airlie, VA. 

7/85, S.T. Wu, To attend AIAA 18th Fluid Dynamics and Plasmadynamics and 
Laser Conference and chair Technical Program Commitrtee, Cinncinnati, OH 

10/85, S. T. Wu to attend CPP Meeting at HAO on Coronal Dynamics in Boulder, 
Co. 

10/85, S. T. Wu to visit Lockheed to discuss Space Lab II results and visit 
China University regarding Solar MHD Simulation Method. 

7/86, S. T. Wu to attend XXVI COSPAR; SMA and SCOSTEP and Meudon Observatory 
8/86, S. Panitchob to visit NOAA/SEL to work on numerical code. 
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